The Random Field Ising Model on Hierarchical Lattices I: 
Phase Diagram and Thermodynamics 



O 

o 

(N 



C/3 



Alexandre Rosas and Sergio Coutinho 

Laboratorio de Fisica Tedrica e Computacional, 

Umversidade Federal de Pemambuco, 

50670-901, Recife, Pemambuco, Brazil. 

The phase diagram and the thermodynamics of the random field Ising model (RFIM) defined on a 
family of diamond hierarchical lattices of arbitrary dimension and scaling factor fe = 2 is investigated. 
The phase diagram is studied considering the flow of the renormalized joint probability distributions 
of couplings and fields. A continuous (Gaussian) and a discrete (delta-bimodal) initial symmetric 
probability distributions for the random fields with variance Hq are particularly considered. The 
thermodynamics properties (energy, specific heat and magnetization) are obtained by an exact 
recurrence procedure and analyzed as function of the temperature and the random field. Close 
and above the paramagnetic-ferromagnetic transition evidences of the formation of rare but large 
field induced correlated clusters of reversed spins are observed. For all studied properties no strong 
qualitative distinct behavior is found whenever the continuous or the discrete distribution of random 
fields are considered. 
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I. INTRODUCTION 

The study of random systems with quenched disor- 
der deserved many experimental and theoretical efforts 
in the last quarter of century, but no significant progress 
has been achieved to full understanding its phase tran- 
sition and critical behavior. In general the bond or the 
site disorder are considered to modeling real systems, the 
paradigms of these models being the Edwards- Anderson 
spin-glass rupdel and random field Ising model (RFIM), 
respectivelytl. Regarding its critical behavior, the RFIM, 
which is believed to belong to the same class of universal- 
ity as the|.diluted antiferromagnet under a uniform mag- 
netic fielda, is known to undergo a phase transition from 
a disordered phase to a long-range ferromagnetic phase 
when the temperature and strength of the random field 
are decreased below some critical values for systems with 
dimension greater than 2. However, the nature of such 
transitions remains not well understood even whether 
it is first-order or continuous. Despite many difficul- 
ties arising from non-equilibrium effects, recent experi- 
ments on magnetically diluted antiferromagnets, such us 
FexZni-xF2, with high concentrations and under high 
fields, provide an experimental characterization of the 
equilibrium critical behavior of the Sd-RFIMB. The lack 
of exact solvable models beyond mean filed-like approx- 
imation has driven the recent theoretical efforts to use 
sophisticated numerical simulation approaches to study 
the temperature and field behavior oLthe. energy, specific 
heat and magnetization of the RFIMqu'EI. The phase di- 
agram has been -also subject to investigation by Monte 
Carlo simulationEl. In this work we investigate the ther- 
modynamic behavior of the RFIM defined on diamond 
hierarchical lattices (DHL's hereafter) with scaling fac- 
tor 6 = 2. 

Hierarchical lattices have been used, since long time 
ago, as a framework to study quenched disordered mag- 



netic systemsQ, as welL-as to study homogeneousQ and 
deterministic aperiodict2l spin models. Whenever viewed 
as a scheme of real space renormalization approxima- 
tion they offer only crude approximations for the critical 
properties of the corresponding homogeneous systems on 
Bravais lattices. However, for random systems like spin 
glasses, they give results comparable to other approaches 
for three-dimensional lattices, such as Monte Carlo sim- 
ulations. It is worth to mention that the values of the 
critical temperature and of the critical exponents asso- 
ciated to the order parameter and correlation length for 
the Ising spin glass model on diamond hierarchical lat- 
tices with graph fractal dimension dp ~ i are surpris- 
ingly close to those obtained by numericaJ— simulations 
for the model defined on the cubic latticdlll. It seems 
that the randomness of couplings and/or fields washes 
out the effects caused by the lack of translational in- 
variance and by the concomitant high inhomogeneous 
coordination-number distribution of the hierarchical lat- 
tices. On the other hand, when viewed as an exactly 
solvable model approach for random spin models, they 
can provide exact results to analyze the dependence of 
the order parameter on temperature and field. For in- 
stance, the local structures of the Edwards- Anderson or- 
der parameter of the short-range Ising spin-glass model 
defined on thcJDHL's has been analyzed by means of this 
methodologyllj. Moreover, concerning the critical behav- 
ior of the latter model, it was possible to obtain. strong ev- 
idences of absence of breaking of universalitytStj regard- 
ing the probability distributions of coupling cpnstants, in 
contrast to what is claimed by some authorsEj. 

Here, we adapt to RFIM the methodology developed 
to the Ising spin-glass model under a fieldEEl. Our aim in 
this paper is to study the phase diagram and the temper- 
ature and field dependence of the thermodynamic poten- 
tials (internal energy, specific heat and magnetization). 
The phase diagram is explored by the renormalization- 



group point of view. The thermodynamic properties are 
obtained from the site-to-site magnetization, which by 
it turn is calculated via a combined method encompass- 
ing real space renormalization-group approach and an 
exact recurrence procedure. We consider both a con- 
tinuous (Gaussian) and a discrete (delta-bimodal) initial 
symmetric probability distributions for the random fields 
with variance Hq. In a subsequent paper, we focus our 
attention to the ground state critical properties by ex- 
tending the present approach to the T = limit. 

The remaining of this paper is organized as follow. In 
the next section we define the model and discuss the 
renormalization flow. In section UIU we analyze the lo- 
cal magnetization either with emphasis in the thermody- 
namics properties or in the local magnetization structure. 
Finally, in section pv|, we summarize our results. 



II. THE MODEL HAMILTONIAN AND THE 
RENORMALIZATION FLOW 

We consider the RFIM on a DHL with an arbitrary 
fractal dimension. This family of lattices is recursively 
constructed starting from a basic unit (called first gen- 
eration) composed by two roots sites connected by a set 
of p parallel branches, each one formed by a series of 
b connected bonds. Further generations are obtained 
by replacing all bonds of the previous one by the ba- 
sic unit itself, as sketched in Figure n^. The resulting 
lattice is a self-similar graph with the fractal dimen- 
sion dp — I + logp/log5. The number of sites and 
the number of bonds are respectively given by Ns = 
2+ {b- l)p[{bp)^ - \\l{bp - 1) and Nb = (bp)^ , N 
being the number of generations or hierarchies. 

The RFIM energy is defined by: 
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where Jij = Jq > are the ferromagnetic coupling con- 
stants of the pairs < i,j > of nearest neighbors spins, 
{Hi] are the quenched local magnetic field at site i and 
(Tj = ±1 are the corresponding Ising spin variables. The 
random field variables {Hi\ are chosen from a probability 
distribution P{H) such that [H'^]h = H^ and [H]h = 0, 
[■■]h meaning the average over the random field disor- 
der. In the present work we consider the Gaussian and 
the delta-bimodal distributions given respectively by: 
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Building a N-generation lattice with uniform ferro- 
magnetic configuration of couplings constants and a ran- 
dom configuration of local fields defined accordingly to 
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FIG. 1: Schematic construction of a diamond hierarchical lat- 
tice with a basic unit composed by p branches and (6 — 1) 
internal sites (scaling factor b). 



a chosen distribution, one can decimate the spins vari- 
ables introduced in the last generation obtaining a (N-1)- 
generation lattice with a new configurations of correlated 
effective couplings and local fields. These new coupling 
and fields are indeed locally correlated accordingly a joint 
probability distribution V{flij), where flij = {Jij, hi, hj} 
represents the set of variables triplets associated with a 
given renormalized branch of a given basic unit of the 
previous generation after the decimation of its (6—1) 
internal sites, as sketched in the Figure 0. 

Nevertheless, to obtain the new random variables of 
the renormalized lattice one has to take into account 
the contributions coming from all decimated branches 
following the lattice topology. Each new coupling con- 
stant being just the sum of the renormalized coupling of 
each branch belonging to the corresponding previous ba- 
sic unit, while the field variable associated with a given 
site being the sum of the renormalized fields coming from 
all branches (of all basic unities) arising from that site, 
in addition to the pre-existing field acting on that site. 
Such procedure weight correctly the sites accordingly to 
the high inhomogeneities found in the site coordination- 
number distribution along the lattice. The decimation 
process of a given basic unit gives rise to the following 



renormalized variables as usual. We have also defined 
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FIG. 2: Renormalization process of a single branch to gener- 
ate the joint probability distribution 'P(Slij). 



renormalization equation for the coupling and fields: 
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where /C = (3 J, hi = fiH and the prime refering to the 
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Here, max(x,y) and min(a;,y) are functions that return 
the maximum or minimum of their arguments respec- 
tively. 

Equations are rather general and were firstly deduced 
in the context of the short range Ising spin glass (SRISG) 
model under an external uniform magnetic fieldllj. Ac- 
tually, for both systems, RFIM and SRISG model, equa- 
tions ^ govern the evolution of the joint probability dis- 
tribution of couplings and fields under the real space dec- 
imation scheme, the distinction between the two models 
being solely introduced by the choice of the initial dis- 
tributions of coupling and fields, which determines the 
evolution towards the corresponding "fixed-point" distri- 
bution. Within the renormalization approach, the phase 
diagram of the RFIM can be directly determined by the 
overall behavior of the flow of the joint distribution un- 
der the renormalization process. To generate this flow we 
adopt the numerical reservoir procedure where pools of 
correlated coupling constants and fields {^ij) are build 
up from the equations ^. For zero fields, this approach 
recovers the one often used to study the flow of the 
coupling constant distribution of the zero field SRISG 
modelQ. In that case, it is worth to mention that the 
distribution evolves within the manifold defined by the 
trivial "zero-field distribution" (a delta-function at the 
origin), within the general space of the joint distribution 
the couplings and fields. For the RFIM case, we consider 
the initial distributions for the fields as given by equa- 
tions H and by a delta function P{Jij) = 6{Jij — J), 
J > 0, for the coupling constants. We analyze the evolu- 
tion of the first moment of the renormalized distributions 
of couplings as a function of the temperature and the 
variance Hq of the initial distribution of random fields. 
The nature of the magnetic phase at a given tempera- 
ture and field Hq is determined by monitoring the evolu- 
tion of J = [J] J and and the field distribution variance 
H = [H^]^/"^, for a iV-generation lattice. As TV ^ 00 the 
paramagnetic phase is stable provide limjv— ►00 J — *■ and 
liniAr^oo H -^ 00, while the ferromagnetic phase is stable 
whenever liniAr^oo J — > 00 and liniAr^oo J/H — -> cx). 

To obtain the phase diagram of the RFIM on a d^ = 3 




FIG. 3: Typical example of a flow diagram for a dp = 3 DHL. 
For a Gaussian distribution with T — 0.8J/fcs, the transition 
field can be determined as the one for which the behavior 
of J/H changes from the paramagnetic to the ferromagnetic 
behavior. The following fields are plotted from the top to the 
bottom of the graphic 0.6, 0.7, 0.8, 0.856, 0.858, 0.9, 1.0, 1.1 
and 1.2. 



lattice we consider a diamond hierarchical lattice with 
6 = 2 and p = 4 and analyze the evolution of the joint 
distribution accordingly to Eq. 0. In figure we show a 
typical plot of the evolution of J/H for a fixed temper- 
ature and several values of Hq up to A'^ = 30 renormal- 
ization steps. The phase boundary can be traced in the 
diagram Hq/Jq against ksT/Jo, as shown in figure ^ for 
the di? = 3 and the dp — 2.58 .. . hierarchical lattices. 
There are slight differences at low temperatures when the 
plots obtained by considering the Gaussian or the delta- 
bimodal distribution of fields are compared, giving rise 
to distinct zero critical fields for the Gaussian {H§ ~ IJ) 
and delta-bimodal {Hq =~ 0.9 J) models. However, both 
plots give the exact zero field critical temperature for the 
pure ferromagnetic Ising model on this lattice given by 
ksTc/J = 1.38... for the three dimensional lattice. In 
Fig y we also plot boundary for the region where the ac- 
tion of the random field is more remarkable. This bound- 



ary is determined and discussed below in section HIE . 

It is worth to mention that we have also performed a 
phase diagram for the dp = "2 DHL, but no ferromag- 
netic phase was observed - even for very small fields, the 
renormalization flow goes to the paramagnetic phase. 



III. THE LOCAL MAGNETIZATION 

The method for evaluation of the local magnetization 
is based on the equivalence of the partition function of 
the model Hamiltonian for a lattice with N generations 
with the one of a single basic unit introduced in the iV*'* 
step with effective fields acting on the spins of its root 
sites in addition to an effective interaction coupling these 
spins. These unknown local effective interaction and 
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FIG. 4: Phase diagram Ho/ Jo against T/Jq, for a d^ = 3, 
open circle (Gaussian) and stars (delta-bimodal), and dp = 
2.58..., open squares (Gaussian) and x (delta-bimodal), 
DHL. The open diamonds (Gaussian) and the crosses (delta- 
bimodal) are an approximate boundary for the region where 
the random field act strongly. 



fields mimic the influence of the remaining lattice spins 
couplings transmitted by the root spins of that unit. Ac- 
tually, these effective interaction and fields arise when we 
perform the trace over the spins variables of the remain- 
ing lattice considering as a decoration unit coupling the 
root sites in the spirit of the decoration transformationliZI. 
The local magnetization and the pair correlation function 
of all sites of the basic unit of the effective system can 
be calculated as a function of its coupling constants and 
local fields and the unknown effective coupling and fields. 
From a subset of these equations, the unknown variables 
can be calculated as a function of the local magnetization 
and the pair correlation function involving the internal 
sites. Now, the unknown effective fields and couplings 
can be eliminated ending up to coupled recursive rela- 
tions giving the local magnetization and the pair corre- 
lation functions involving the internal sites in terms of 
the corresponding values involving the root sites. This 
approach has been extensively applied to investigate the 
local order parameter of several magnetic models defined 
on hierarchical lattices. Those studiei s i co mprise the fer- 
romagnetic rUing model with uniformllSilH jand aperiodic 
interactionscJ, the-spin-glass Ising modeltda and the fer- 
romagnetic PottsEil model defined on diamond hierarchi- 
cal lattices, as well as the Ising model (both the pure and 
the spin-glass cases) defined on the Wheatstone bridge hi- 
erarchical latticeEj and on the m-sheet Sierpinskii Gasket 
fractal latticeEj. 

For the RFIM the coupled recursive equations between 



local magnetization and pair correlation functions are 
written as 



< a > = Ai + A2 < ^ii > 

+ Ag < fi2 > +^4 < ^J■lfJ.2 > (4a) 

< cr/ii > = A2 + Ai < ^i > 

+ A4 < fl2 > +A3 < ^1^2 > (4b) 

< cr/i2 > = A3 + ^4 < Ml > 

+ Ai< fi2 > +A2 < Aii/i2 > (4c) 

where a is the magnetization of the internal site of a given 
branch, /xi and /i2 are the magnetizations of the root 
sites of this branch, h, Ki and K2 and the interaction 
constants acting on a. The coefficients of equations (0) 
are given by 



AAi = tanh(/i + Ki + K2) + tanh(/i + Ki- K2) 

+ tanh(/i - Ki+ K2) + tanh(/i - Ki~ K2), (5a) 
4yl2 = tanh(/i + JsTi + K2) + tanh(/i + iCi - ^^2) 

- tanh(/i - iCi + is:2) - tanh(/i - isTi - 1^:2), (5b) 
4A3 == tanh(/i + JsTi + ^5:2) - tanh(/i + iCi - ^^2) 

+ tanh(/i - Ki+ K2) - tanh(/i - isTi - K2), (5c) 
4A4 = tanh(/i + Ki+ K2) - tanh(/i + iCi - K2) 

- tanh(/i - iCi + ^5:2) + tanh(/i - A'l - K2). (5d) 

For zero field systems, A3 = A4 = and equa- 
tions decouple. In this case, the equation Ea was 
used to study the critical and multifractal properties of 
the Edijpj-ds -Anderson order parameter of the SRISG 
mode]|13o and t]ie,|ttjagnetization of the pure ferromag- 
netic Ising modeHEJ 

To obtain the site-to-sitc magnetization and the bond- 
to-bond correlation functions of the whole lattice we pro- 
ceed the following script: firstly, we consider a N gener- 
ation lattice and successively decimate the spins up to 
the first generation (basic unit), keeping track the cor- 
responding local values of couplings and fields variables 
(fiij) at each step. Then, we calculate the initial values 
for the magnetization and pair correlation function of the 
root sites, using the effective partition function. Now, we 
recursively calculate the local magnetization < ct^ > and 
the bond correlation function < ui^, > using at each step 
the exactly stored local values for Q.ij . This procedure re- 
quires large computational efforts as the number of sites 
and bonds grow exponentially of with the number of gen- 
erations Ns ^ Nb ^ h^ . Moreover, one has to average 
over many realizations of the initial distributions P{H) 
(samples) , yielding to the magnetization per spin and the 
dimensionless internal energy per spin straightforwardly 
obtained by 




FIG. 5: Average magnetization as a function of the tempera- 
ture and the strength of the random field for the RFIM with 
a delta bimodal distribution and dp ~ 3. 
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where [...]c denotes the configurational average taken over 
many samples and hi = Hi/Jo- 



A. Thermodynamic properties 

In figures ra to 0, we present the temperature and field 
dependence of the magnetization, the internal energy and 
the specific heat, obtained from equations ra and m - the 
specific heat being calculated by the numerical derivative 
of the internal energy with respect to the temperature. In 
each figure, we show the results for either the Gaussian or 
the delta-bimodal probability distributions calculated by 
the exact methodology considering a d^ = 3 DHL with 
N — 7 hierarchies which correspond to Ns ^ 1.2e + 6 
and Nb ^ 2.1e -I- 6 sites and bonds respectively. 

In Fig. |5| and |6|, whatever is the probability distribu- 
tion, the magnetization decreases faster for high fields 
(below the critical fields) as the temperature is increased, 
indicating that stronger is the field greater is the prob- 
ability of appearance of correlated clusters of reversed 
spins with respect to the magnetization. A detailed anal- 
ysis of the distribution of reversed clusters, at zero tem- 
perature, as a function of the random field and an ex- 
ternal uniform field is now under consideration and the 
results will be present elsewhere. Figures ^ and g also 
show that the field and temperature dependence are qual- 
itatively similar with respect to the random field initial 
distribution. In the subsequent paper, we analyze the 
dramatic behavior of the magnetization at T = 0, sug- 
gesting a discontinuity with /S — exponent. 

In figure Fig. [^, the behavior of the internal energy 
as function of the temperature and field strength is pre- 
sented for the Gaussian model (the delta-bimodal one 
being very similar). The internal energy increases with 
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FIG. 6: Average magnetization as a function of the temper- 
ature and the strength of the random field for the RFIM with 
a Gaussian distribution dp =2). 



FIG. 8: Specific heat as a function of the temperature and 
the strength of the random field for the RFIM with a delta- 
bimodal distribution and dp = Z- 
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FIG. 7: Internal energy as a function of the temperature and 
the strength of the random field for the RFIM with a Gaussian 
distribution and dp ~ S. 



FIG. 9: Specific heat as a function of the temperature and the 
strength of the random field for the RFIM with a Gaussian 
distribution and dp ~ 3. 



the increasing of the temperature for a fixed value of the 
field the slope becoming bigger (smaller) for values of the 
field below (above) H§{T). For small temperatures and 
fields, the internal energy is dominated by the exchange 
term, so that it is almost constant in this region. For high 
fields, however, it is the local field term that dominates it, 
thus justifying the linear dependence the internal energy 
upon Hq, for Hq:^ H§. 

The specific heat behavior is displayed in Fig. g for the 
delta-bimodal model and in Fig. O for the Gaussian one. 
We notice that the temperature where the maximum oc- 
curs almost does not depend upon Hq, and consequently 
it does not coincide with the corresponding transition 
temperature Tc{Hq). This behavior is characteristic of 
Ising spin systems defined on t hs d iamond family of hi- 
erarchical lattices with dp > 2li3o. For fields near and 
above the critical one, -ffo(O), we observe that the specific 
heat for the delta-bimodal model displays a "plateau". 



absent in the Gaussian model. This difference appears 
only in the specific heat since, among the studied ther- 
modynamic potential, it is the most susceptible to fluc- 
tuations. 



B. Local magnetization structure 

In this section, we focus our attention on the struc- 
ture of the local magnetization, in order to analyze the 
effect of spin reversion. The Edwards- Anderson parame- 
ter qea = l/J Vgy^ „-[< (Ti >^]c is calculated and shown in 
figures |l^ and |ll|, for the delta-bimodal and the Gaussian 
model respectively. In both cases, there is a depression 
the transition for low temperatures and field near and 
above the critical one. This depression is a consequence 
of the competition between the exchange and the field 
terms, which in the transitions are almost equal, giving 
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FIG. 10: Edwards-Anderson parameter as a function of the 
temperature and the strength of the random field for the 
RFIM with a delta-bimodal distribution and dp ~ 3. 



FIG. 12: Variance of magnetization distribution as a function 
of the temperature and the strength of tire random field for 
the RFIM with a delta-bimodal distribution and dp = 3 . 
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FIG. 11: Ed wards- Anderson parameter as a function of the 
temperature and the strength of the random field for the 
RFIM with a delta-bimodal distribution and dp = 3. 



rise to small local magnetization. It is worth to notice 
that the depression is more evident in the delta-bimodal 
model, specially for low temperatures, since the probabil- 
ity of exact frustrated spins are non zero in this model. 
For higher fields, the spins follow their local fields and 
consequently qea —^ 1- Analogously, for small fields and 
temperatures (ferromagnetic phase), all spins are corre- 
lated and aligned, so that we also have qea — > 1- 

In order to investigate the role of the reversed clusters 
we analyze the temperature and field dependence of the 
configurational average of the magnetization distribution 
Pi{m) (the l-index labeling the samples) variance, that 



A 



dmra Pi{m)\ — [ / dramPi{m)\ 



(8) 



A gives information about the average numbers of site 
with reversed magnetization induced by the random field. 



In figures |ij and |T^ we plot A as function of the temper- 
ature and field, for the Gaussian and the delta-bimodal 
models, respectively. For a fixed temperature, the num- 
ber of sites with reversed magnetization increases with 
Hq, the beginning of the increasing region being deter- 
mined by the corresponding critical field Hq{T). Look- 
ing at figures |l^ and |l^ for fixed fields instead, we notice 
that A displays a maximum above which the temper- 
ature dominates the system behavior and all the local 
magnetizations vanishes. We also observe that the max- 
ima shifts for low temperatures as Hq increases. This is 
a consequence of the decreasing of the critical tempera- 
tures with the field strength. Since the maximum of A 
determines the point where the effect of the random field 
is the strongest, we determine the boundary of the region 
where the field "plays its role" as the point that A is half 
of its maximum. This boundary should be related, with 
the experimentally observed non-equilibrium lineEZI. We 
also point out that for very low temperatures and values 
of Hq close and above Hq one observes a small increase of 
the A parameter, indicating the persistence of sites with 
reversed magnetization at ground sate (this behavior is 
more evident for the Gaussian model). 



IV. CONCLUSION 

We studied the RFIM on DHL with fractal dimension 
three, considering either a Gaussian or a delta-bimodal 
probability distribution for the random fields. We ana- 
lyzed the renormalization flow as well as thermodynamic 
properties. We obtained the phase diagram for the three 
dimensional and for the dp — 2.58 . . . systems. We also 
obtained a boundary line separating the region domi- 
nated by the random fields from the one dominated by 
the temperature, for the three dimensional system. Con- 
cerning to the thermodynamics, we observed that the 
thermodynamic potentials obtained within both proba- 
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bility distributions displays a very similar behavior, the 
Edwards- Anderson parameter at low temperatures (as a 
consequence of exact frustrated sites) and the heat ca- 
pacity near the critical field (by virtue of strong fluctu- 
ations) being the exceptions. The parameter A, which 
gives information about the clusters of reversed spins, 
gives useful information of the random field effect, dis- 
playing useful information about the competition among 
the random field, temperature and the ferromagnetic or- 
der. 
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